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§ 1. Introduction. 

In his " Adams' Prize Essay " # for the year 1856, Maxwell showed that the rings 
of the planet Saturn could only be stable for small disturbances on the theory that 
they were composed of meteorites sufficiently small. This has been confirmed since 
by spectroscopic evidence and is now generally accepted. In continuance of the 
same idea, the various divisions of the rings have been accounted for by presuming 
that, in those positions where a single particle moving in a circular orbit about the 
planet would have a period simply commensurate with that of one of the nearer 
satellites of Saturn, instability would result. This idea has been fully emphasized 
recently by Lowell, f His observations at Flagstaff have disclosed a large number 
of additional divisions in the rings (see Appendix to this paper). They have the 
appearance of fine lines traced on the surface of the rings. In each case Lowell is 
able to show that the divisions occur at intervals of periods commensurable with that 
of satellite Mimas. The periods have the ratios such as f, x T , |r, y\-, &c. Lowell 
has stated the argument for this view in ' Bulletin/ 32, p. 189. If the action of one 
body upon another revolving about a third be examined by the method of the 
variation of arbitrary constants, in the expressions for the periodic inequalities 
in the radius vector and the longitude, there appear terms of the type 
[Q/(pn — qn f )~] cos {{pn — qn f ) £ + Q}, where n and n! are the mean motions of the 
perturbing and perturbed bodies, p and q are integers, and the remaining quantities 
are constants. It is clear that when the ratio njn f is approximately equal to q/p, 
then the inequality will become very large. 

We may take a satellite of Saturn as one of the bodies and one of the particles 
forming the ring as the other ; if n\n! == q/p, approximately, then the particle will 

* Maxwell's ' Collected Works/ I., p. 288. 
t Lowell, ' Observatory Bulletin,' No. 66. 
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depart considerably from its unperturbed path and collision with other particles will 
result. In this way the divisions in the ring have been explained. 

Some doubt has been cast upon this theory, and it has been shown # that even 
when n and n r are commensurable, a closer examination of the motion leads to the 
conclusion that the denominator will not vanish. 

It is also noticeable that this explanation takes no account of the attraction of the 
numerous particles upon one another, which may be considerable. 

A re-examination of the matter is made in the present paper. As the satellites of 
Saturn are all approximately in the same plane as the ring, the problem is formulated 
in two dimensions only. The satellite is assumed to follow an unperturbed circular 
orbit, and the problem reduces to a slight variation of the " restricted problem " of 
three bodies. We shall consider the effect of this satellite upon a number of particles 
forming a single ring round the planet, subject to their mutual attraction as well as 
that of the satellite and of Saturn. The actual Saturnian rings are supposed to be 
composed of a number of such rings arranged concentrically. These will have some 
effect one upon the other, but, for the present, this effect is disregarded. 

In his paper, Maxwell considered the single ring of particles only. He found 
that the equations of motion could be satisfied by assuming that the particles rotated 
round the primary in a circle with suitable angular motion. He then examined the 
effect of a small arbitrary disturbance upon them, and his results show that the 
disturbances would remain small if the masses of the particles were sufficiently small. 
That is, the ring would be " ordinarily" stable. 

In the present paper the plan is different. The disturbance of the ring of particles 
by the satellite is examined, with a view to determining under what conditions the 
departure from a certain fixed circle will be large. It is clear that if the departures 
do become large, collisions with adjacent rings of particles will result, and the particles 
will leave the vicinity of the original circle irrevocably. In this case a division in 
the ring will result. It is with this meaning that the terms stability and instability 
have been used in the paper. But, as will be pointed out again in its proper place, 
the orbits in which the departure from the circular form does not become great with 
increase of time may yet become " ordinarily " unstable if further small arbitrary 
displacements are imposed upon them. 

The results of this paper will therefore indicate some, but not necessarily all, of the 
positions of divisions in the rings due to instability of whatever kind. 

In §§ 2 to 4 an analytical theory is fully worked out on the supposition of equal 
particles in each ring. In § 5 it is shown how amendments may be introduced to 
cover the case of unequal particles. The application to the Saturnian system is given 
in § 6, and the last paragraph summarises the results obtained. 



* Tisserand, ' Mec. Celeste,' vol. iv., p. 420. 
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§ 2. Formation of the Equations. 

Let M be the mass of the primary and m! the mass of the principal satellite which 
is assumed to describe an unperturbed circle round the primary. Take the origin at 
M. Let there be n particles forming a ring round the primary, subject to attraction 
from M, m\ and one another, and let the mass and co-ordinates at time t of particle X 
be m x , r A , (9 A . If the co-ordinates of m f at the same time are r\ 6\ then the motion 
of particle X will be produced by forces which are the derivatives of the function 
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As we are assuming that m! describes an unperturbed circle, 



r f = a' 



and 



(9 = CO t~\-€ , 



where W 3 = M + m' = M, with sufficient approximation. 

Let us assume now that the remaining particles are moving in the vicinity of the 
vertices of a regular polygon of radius a. Then we may put 

r A = a-f-px, 

A = otf-fe + A . 2tt fn+.a-\, 



for all values of X from 1 to n, where p and <r are assumed small, so that squares, 
products, and higher powers of them and their first derivatives with regard to the 
time may be neglected. 

The equations (l) now reduce to 
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In order to determine the derivatives, we write in the formula for A A; ^ A -- f ~6 x , 
and a = r A /V. Then 

A A ~ X = {l+a 2 -2acos</> A }~*-rV 

= {|f& -1- &i cos ^ + ... + 6* cos ^ + . . . } -f- r' 

by Fourier's series. 

This series will be taken as absolutely and uniformly convergent. 
We find then 



3P 

3r A 



M + m A m! \ t db 



2 



+ 



72)2 



3 



+ . . . 4- 



OL 



w cos %<h + . . . r ^ c os 

da r 2 



2m 



/* 



M 



V A -7ycos((9 / ,-0 A ) cos(0„- 

IV V 



0.) 



3 2 F 



'/V 



/x 3r A 3r 



H- 



3 2 F 



/0-, 



% "3ft,3r, 



fJL^ J ' A 



^3F 

n 3#a 



m 



n M+m A , m 3 W /l7 , 7 • \ 

2 5 h — —-2 (f 6 + ... +6.- COS £0 . . . ) 

r A 3 r' 3 3a ^ * r 7 



2m, 



'* \J^ 3 






/5 A 



+ 2m 



^ 



/* 



co^(^-A) + 3_{r k -r^^ 



Ukfx' 



T) 6 



+ 



2cosj^-e A y 



^ 



^j 



m/ 3 



772 



r /2 3a 

— 2m 



/7 * *7 ' * \ lit • 

{o l sm ^ + ... + *6,- sin *0 + . . . ) tj sm 96 



r 



o*x 



^ 



r^ sin (0 



/x 



0x) __ 3 {r A -r M cos (d^ —dx)} n r, sin ( ^--(9 A ) 



I) 3 



T) 5 



_ sin _(e^ L - 6x) 



p 



(cr^ — orx), 



m 



It 7 -7 \ 7/t • x , 

j— { . , , + *o t - sin *0 + . . . j ^ sm + 2m^ 



ry+' ry+ 



? ysin(g /x -0 A ) 



D 3 



-2 sin (O^-Oa) 



/^ 



SATELLITES UPON THE FORM OP SATURN'S RING. 



105 



3 2 F 



^ p " ndr.de, 



ml 3 



r'V x 3 a 



(b l sin <£+ ... -\-ib i sin i(p+ ...) 



-. — - 9 ( . . . + ib { sin id)...) 

r f r* x 



— 2m 



H- 



h>- 



f>\ 



-t-Zm 



p 



^ 



3ry sin (Op- Ox) {JjlZLl cos ( /f-0x)} 

sin^-flx) _ Sty sin^ -A) {r ^-r A cos (^-0 A )} « 

,3 " " ~™"~ f) 5 



D; 3 



H 3 sin (0 



p> 



$>) 



& 



v» 



3 2 F 



* r A 3^30 A 



m ; / . 27 . \ , m/ 

—r— ( . . . r&i cos ^0 . . . ) + — t^ cos <h 

r'r k x T J r f2 T 



<r\ 



+ 2m 



& 



h- 



V ft cos(^- -e A ) 3r/r A sin 2 (^-0 A ) 



-*~^/xA 



T) 5 



-COs(^-# A ) ((T^-o-x). 

In the summations of the right-hand members, /x takes all integral values from 
1 to n 9 except jm = X. 

The zero values of these derivatives are obtained by putting 
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In the summations of the right-hand members, jul takes all integral values from 
1 to n, except jul = X. 

Next assume that all the small particles forming the ring are equal to one another. 
That is m A = m. 
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The quantities K, L„ M„ ft, can readily be found by direct summation when n, the 
number of particles, and s are known. 

In re-writing the differential equations (2), we may now omit the suffixes of p and <x. 
Change the independent variable from t to <f> = (<■/-«) £ + e'— e— X . 2?r/n. Also put 
m/M = v, m ; /M = i/, o)7« = *, (ic , -l)- 1 = ic, and, to secure homogeneity, replace 
P by ap. Let us further assume that co 2 a 3 = M, and «'W 3 - M (the latter holds very 
approximately when ra' describes a circle), so that we have a/a' = (w/w) 2! \ 

The differential equations then become 
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The equations (2) may be replaced, under the suppositions made, by equations (3). 
Equations (2) form a system of n pairs of linear equations of the second order. The 
complete integral will therefore involve in arbitrary constants. The system of 
equations (3) will give the same result, for the solution of (3) will be a function of s 
involving four arbitrary constants. By giving s its n values, 0, I, 2, ..., (n— l), we 
arrive at the complete integral involving in arbitrary constants. 

Now it has been shown by Tisseband # that for large values of n, whatever the 
value of s may be, the limiting value of L s is 0'0194^ 3 , N s = 2L S and M s = 0. 

These values largely simplify the discussion of the stability of the system. 
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'M^c. Celeste/ vol. ii., p. 184. 
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The best methods of determining the values of 6 r , and its first and second 
derivatives for known values of a, or a/ a', are given by Tisserand. # The complete 
evaluation of a number of these quantities for various ratios, applicable to the solar 
system, is given by Ponteooulant.! For the purpose of estimating the order of the 
numerical values of the quantities r>s , we may take the highest ratio a likely to 
occur as that of the outer edge of the ring to the mean distance of Mimas. This 
ratio is 07461 (see Appendix for data). Pontecoulant gives the values for 
a = 072333, which we may use to avoid laborious calculation. If we take / = 7 . 10~ 8 , 
the value for Mimas, we find 

9 lf0 = -20*2590 [6'7530vL s ; 

9 lfl = -1'3'5. 10-°, O h2 = -l"53..10- 6 , e if 3 = -l'63-. 10- 6 ...; 

9 2>1 =--4'07 . 10~ 7 , 2 , 3 = -7'39.1O- 7 , 2>a : = -9*34 . 1CT 6 ... ; 

3>o = -1-46 . 10" 7 -179 . 10- 7 K, 9 3fl - 4*07 . 10" 7 , 9 3 , 2 - 3*69 . 10" 7 .;. 
e 5 , = -13'5060i/L,. 

It is clear that, compared with O ]?0 , all products and squares of the remaining 0's 
may be neglected. 

§ 3. Solution of the Equations. 
(a) The complementary function. 



The equations 

p n — 2/ccr / + p20 1)r cos r0 + <r29 2ir sin r<j> = 0, 
<r f + 2/c/) / + / o29 4)r sin ?v/> + o-26 5) r cos r<p = 0. 



(6) 



belong to the class of homogeneous linear differential equations with periodic 
coefficients. The integral is kiiown to be the sum of the forms e c(f> f (</>), where /(</>) 
is a periodic function of <j> with the same period as the coefficients in the equations (6). 
Equations of this form in one dependent variable have been discussed by Whittaker,J 
Young,§ Inoe,|| and Baker. ^ The present solution is a simple extension of the work 
of these writers. 

Let 

P = e cr/, A, 

a- = e c< % 

* ' Mec. Celeste,' vol. i., p. 270, et seq. 
t < Systeme du Monde,' vol. 3, pp. 353-376. 

X 'Proc Inter. Congress Math.' vol. 1, 1912; ' Proe. Edin. Math. Soc.' xxxii., p. 76. 
§ * Proc. Edin. Math. Soc' xxxii., p. 81. 
|| 'Monthly Notices E.A.S.' lxxv., 5, p. 436. 
1[ H. F. Baker, 'Phil. Trans.' A., vol. 216, p. 129. 
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where A and X are, as has been* said, purely periodic functions of period 2w. On 
substituting in equations (6) we find 
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On eliminating A and X we find 
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Following the method of Whittaker previously referred to, let us now assume 
that 
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or 

A 2 2 

95,0 rh 

* The use of n is to be distinguished from a former use where it referred to the number of particles in 
the ring. 



SATELLITES UPON THE FORM OF 'SATURN'S RING. 



Ill 



We now substitute the assumed values for A, X, c, ljO in equations (7) and equate 
to zero the coefficients of each term in 6 r> „ ?v? , Ptq9 &c. It will be found that the 
relations (8) are satisfied identically. Two conditions further must be imposed in 
order that all the unknown coefficients may be determined. These are : 

(i) The term cos (n<p~r) must not appear in the series for A ; 

(ii) The solutions for A and X must be purely periodic with period 2ir. 

The condition (i.) amounts to a definition of r, and condition (ii) secures that no 
part of the exponent shall appear in the periodic series. Further, these conditions 
determine uniquely the undetermined coefficients in the series for 1>o and c. The 
work from this point is purely mechanical though long. The following sample 
sufficiently indicates its character. 

On equating to zero the terms involving ]>r we find 

- 2c h r nA cos (n<j)—r) -f A" h r — 2kc u r X cos (n<j>—r) 

— 2icX f hr + a 1A) A hr + a hr A sin (n^ — r) + A cos r<j> sin (rup — r) = f (12) 
— 2c lt r nXo sin (n</> — r) + X" lt r + 2icc h r A sin (n<f> — r) + 2k A' h r + 5j X 1? r = 



In the case when r is not 2w or n, it is clear that 



C l,r 







and 



a 



1, r 



0. 



Equation (12) then reduces to 

A\ r — 2icX , 1>r H-a lt0 A ltr + |-A [sm {(n-f r) 0— t} + sin {(n—r) <£~-r}J = 
X 1 r -f- 2/c A j 5 r + 5j X 



-1, r 



0. 



► • ( 13 ) 



Solving in the usual way we find 



.-A-l.r 



X 1>r 



A n 2 {{n + r) 2 -O^ } sin {(n + r) < f>-T} A 7i 2 {{n-r Y-er^} sin {(n-r) 4>~r] 
2r(2n + r) (a lfO 5tO — n 2 {n + r) 2 ) 2r(2n—r) {a^ Q^ — n 2 (n-r) 2 } 

A u 2 k (n + r) cos {(n + r) ch—r) . A u 2 k (n—r)- cob {(n—r) (j>—t} 

f . . — - ? ~ — ; — ■ — — t — ■ — J— — < ' ' ■ ' ' ; • 

r{2n + r) {ci X0 Q b , -~n 2 (n-i-r) 2 } r(2n—r){a h0 Q 5}0 —n 2 (n—r) 2 } 



In the special case where r = n, we have 



and 



'1,7? 



o, 



<h,n = 0. 



A /; liW — 2ic,X' 1>w -l-a lt oA lt7? + i-Ao {sin (2n<p — r) — sin r} 
X' v>1 -f 2/c A x 7? + 5j . X ls v 



0. 



"^ 



> . . . (14) 



From which 



and 



A = Ao(4w 3 -9 5 , )sin(2n 9 !,-T) A . 

6(4w 4 -a li0 e s ,o) 2a lj0 



X, 



n 



2 /C^Ap COS \2U(J) — t) 

* 4w 4 -a lf0 e 5f0 

B 2 
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Again, In the special case where r = 2^,, we find in place of equations (12), the 



following 



2c h2n nA cos (n<f> — r) ~\~ A" h2n — 2KC h2n X cos.(n<f> — T) — 2/cX / 1 , 2n 

+ a lj0 A liSn 4'« li2 »-A^ sin (rujy — r) -!•-§- A {sin (3^--~t) — sin (nc/i — r) cos 2r 



^ 



cos (uip — r) sin 2t; 



\ — 



0. 



S> (15) 



■ — 2c lf 2w nX sin (nc/> — t) -f X" lt 2w + 2/cc 1? Sb A sin {71$ — r) ~h 2/cA/ ]? 2n + 6t X l5 2n , = 0. J - 
In order to avoid the explicit appearance of <p, we must have 

Since we have already stipulated that A must not contain any term in cos (ntp — r), 
c h2n must he so chosen as to make quantities involving cos (n^ — r)' annul Hence we 
must have 

{2c li2n nA --2«:c lt2n X --|-A. sin 2r} cos (w,0— t) — 2« / X / 1) j ftt =. 0, 

-2c 1?2? //iX + 2/cc l52w A } siu(^0—T) + X" 1 , 2w -f9 5t0 X 1 , 2n = 0, 



> 



(16) 



Whence 



C l, 2?7 r 



n (9 5?0 — ^ 2 ) sin 2 



T 



4(«i >o 05.o"-^ 4 ) 



v = l. ( n CT i,o) (Qs,o + ^ ) sm 2 tA smj^-T) 
, ^(^,065,0-^) 

To the value for X h2n must be added the further particular solution arising from 
the term |A sin {3n<j>— r) in (15). It is 



1, 2« 



"Y — 

■^-1,292 



(9 5 > o — 9t& 2 ) A q sin (3n^> 
16(a 1>0 e 5(0 -9^ 4 ) 

8 n/c A ,cos (Zn<p—T) ■ 

8(« 1 ,o©rN0-9^ 4 ) 



~-t) 



Proceeding in this way, we have the following results 

Terms not involving argument : 

n A 

A sin (ncp — r). 

X COS (riep — r). 

None. 
2fcnX = ■— (^1,0— ^ 2 ) A , 
a lt0 = n 2 + 4:K 2 n 2 /(Q 6t0 --ri 2 ). 



In X 

JL J- 1. O 

Also 
when 
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Terms involving argument O hr} where r is not n nor 2n : 

I 

» 

A ^ ~ A o ^ 2 {(^ + r) 2 '-9 5 , o }sin{(^ + r)0-T} A ^ 2 {(tt-r) 8 -e 6t0 } sin {(n-r) 0-t} 
1,r "" 2r(2n + r) K 9 5 , ~^ 2 (^ + r) 2 ) 2r (2ri-r) {^, o e 5 , -^ 2 (^-r) 2 } 

_ A n\ (n + r ) cos {(rH-~r) <f>- T } A rA (rc,-r) cos {(^~-r) </ >~-t} 
x ' f + r (2w + r); {a lt e 5 , o-^ 2 {n + r) 2 } ' r (2n-r) {a h ,Q^~~ri 2 (n-r) 2 } ' 

Terms involving argument Q hv : 

A - Ao_(4^ 2 _---e5,o) sin (2u(P-~t) A sin r 
A, '*-~ "6(4^ 4 -a lj0 e 5 , ) 2a li0 ' ; 

— 2/mA cos (2n< j> — T) 
^~~ 3(4^-^,005,0)""' 



Terms involving argument lf2 „.* 

a _ ( °5 p - gjl!) A sin (3n<p - r) 

l>8n """""~"i6K7e fif0 -9w 4 ) 



n 



( 5f p ~-n 2 ) sin 2r 
4KA,o-^ 4 ) 



-/i 



3wkA o c6 s(3 w0-t) (w 2 -«, q) (9,. n -n 2 ) sin 2tA sin jn ^r) 
] ' ? "~" ' 8(ai, e B , _ „-9n 4 ) "' 8Kn(a,. o e B ,o-9/i*) 



Terms involving argument 9 3 , r , w/iere r is not n nor In: 

_ n « {e 5 o _ (n+r) 8 } Xosjn {(w+r)^-T} w 8 {9 6 , - (wh')'}. x p sin {(w..-r)#-£L 

A ^ ~ 'ir{n+r)Jal Q^ -n^n'+r) a } 4r (w-r) {a 1>0 e Bi0 -« 3 (n -r) 2 } 

• X rAc ( w + r) cos { (n±r) <p-r} \ ^ K .h-rl G 9^iil}- r l±I^i 

2 ' r ~ 2r (b+i : )WV : » 2 (« + r) 2 } 2r (»-r) {^.Ar"' (n-r) 3 } 

Terms involving argument 6 2 ,„ : 

A _5y05,o-^ 2 )sin(2^-r)_X 2 _sin2T L 

v- _ _ 2X kw cos (2n$~ t) 
**.«-■- 3(^,095,0-4^) 
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Terms involving argument 



2, 2m 



ku cos 2t /rn 2 sin 2t 



^2, 2» ""* o 2 ? ^2, 2?? ' 



s ,a-» 8 ' ^ 2(n«-a 1 , e 6 . )* 



a ' — (Qfi.p — S^ 2 ) X sin (3n^> — t) 

9w 4 -« 1 ,„0 B>o 



Y _ 3«-raX co s (3n<p— -r) (Q, 4-n 2 ) /mX sin 2t sin {n< p—T) 



Terms involving argument 0^ r , where r is not n nor 2n: 
A _ A n 2 /e (n -f r ) sin { (n + r ) <f> — r} A n 2 /c (n — r) sin {(n—r)<p 



r (2n + r) {a l5O 5 , o -n 2 (n-f-r) 2 } r (2n~~r) {a lf0 ©5,o-™ 2 (^-r) 2 } ' 

X — A n 2 {(n-hr) 2 --a 1<0 } cos {(n + r) <p — r} x A Q n 2 {(n—r) 2 --a M) } cos {(n— r) cp — r} 
4,r ~~ 2r (2n + r) {a h0 6& t0 — n a (n + r) 2 } ' 2r (2n—r) {a X0 Q 5f0 ~~n 2 (n~~r) 2 } 

Terms involving argument 4> „. : 

A - 2n-/cA _ sjn (2^^) — t) 
A ^ 3(a 1 , e s>0 -'4i i r' 

Y — A.'(.4^ il T ty i ") cos (% n 4>_~' T ) _ Ao cos T 
" 4 ' ,! "~ "" 6(a Ii0 e 6>0 -ln i ) ~~ " _ 2e7r- 

Terms involving argument 42w : 

^/c cos 2r n 2 /c sin 2r 



Cvj <>,, — ^ ,-, .« c. 



4 ' 2M e 5 , ~n 2 ' ^ 2B "2(a ]>0 e 6 , - w *)' 



A — 3/c^A sin (3?i(p—r) 

Y — (9^ 2 — ^i o) A cos (3n^>— t) A cos 2t cos (n^—r) 

4,? " " Te K aTo-9» 4 ) 2 (e M -n 8 ) 

(^i, + n 2 ) sin 2tA sin (mp — r) 
4(ai fO 06 f o-^ 4 ) 

Terms involving argument 5> r ^ ?.#/i6n r &s no/: ^ nor 2n : 

A — X n 2 /r (n-f-r) sin {(n + r) 0— t} __ X n 2 yc (n — r) sin {(n— r) ^—t} 
r (2n + r) {a h0 Q 5)0 —?i 2 (n + r) 2 } r (2n—r) {a- lt ' B 5>0 — n 2 (n—r) 2 }' 



v __ X n 2 {(n-ff) 2 - a lift } cos {(n + r) <£— t} X n 2 {(n— r) 2 — a 2 } cos {(n—r) d — 
2r (2n + r) {a h0 Q 5tQ -n 2 (n + r) 2 } 2r (2n-r) {a h0 G- Ot0 -n 2 (n-r) 2 } 



T 



■} 
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Terms involving argument & 5>w : 



^5, n 



o. 



Cb.n — ^j 



A, 



2X,yc2^ sin (2n0 — t) X cos t 



a, «. 



3(a li0 e Bi0 -4n 4 ) 



29 



5,0 



Y — ^o (4n 2 --c& K0 ) cos (2nft — t) 

6(ai, o 06,6--4n 4 ) 



Terms involving argument 5f 2/t : 



^o, 2» 



(^ ~~^ 2 ) cos 2t 
2(B 5 , -n 3 ) ' 



c 



5, 2n 



n-(# lt0 --n 2 ) sin 2t 

4(a 1>0 e 5s0 -n 4 ) 



A ~ 3n/cX sin (3n0-~T ) 
A5 ' 2H ~ 8(9n*-a 1>o e., ) ' 



A 5, 2>i 



9(n 2 —a, t0 ) Xp cos (3n0 — t) (a 1)0 — n 2 ) cos 2tA cos,(^ft— t) 



16(9^ 4 -~6H, o 0o,o) 



■4 



4^(e 5i0 — n 2 ) 



(a h() + n 2 ) (a^o—n 2 ) sin 2tA sin (n<f}—r) 

8/c7^(a uo e 5)U -7^ 4 ) 



Terms involving powers products of the 0's follow in similar fashion. 
If we summarize the parts specially required, we find 

(e M -^ 2 ) (e 5 , -O = 4icV + i(0 5f o-^ 2 ) cos 2T0 ]f2tt --/cn cos 2r0 2 , 2/r f- /m cos 2T0 4 , 2ri 

-i(^i,o-^ 2 )cos2T0 5t2 „+...; . . (17) 



and 

2c (a M O 6 , -^ 4 ) = 2 n (Q5,o-^ 2 ) sin 2T0 lt2w -^n 2 sin 2rQ 2f2u + K?i 2 sin 2t9 4(2?1 

- h n ( a hQ- n2 ) ^ n 2r ©5, 2n + "...;' . 
where, as already stated, 

Ko-w a )(e M -w 2 ) - 4kV. 



(18) 



It is necessary to examine the expressions just obtained in order to see whether 
the complete integral of equations (6) has been found. 

The integer n is determined so as most nearly to satisfy the relation 



when 



(e li0 -n s )(e M -n a ) = 4^ 8 , 

©i,o ? ©5,o an( l K are known. 



The negative value of n will also satisfy this relation. 

On solving equation (17), for each value of n there will be, in general, two values 
of 2t, equal and opposite in sign. So that altogether there are four distinct values 
of 2t obtainable. Each of these with the corresponding value of n will give a 
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different value of c on substituting in equation (18), and different values for A and X. 
Hence there are four distinct solutions and these when multiplied by arbitrary 
constants will give the complete primitive of equations (6). 

(6) The Particular Integral. 

We have now to determine the particular integral of equations (4). We shall assume 
only one general term on the right-hand side and take the complete solution as the 
sum of a series of the corresponding solutions. The equations may therefore be 
written 

p" — 2K<r' + ptQ hr COS r0 + flr20 air sill Tcj> = ^Q^ u fi +in ^ : 

cr"~\~2Kp~\~p^Q A/r sin r<f>-\-cr%Q 6tr cos r<p = 0. 



> . . . . (19) 



Assume 
and 



p = ^Kf 



__ ,,<-ni 
G' — 6 



*x 



where X and A as before are functions of <f>. On substituting in equations (1.9) and 
reducing, we find 



w 



- m 2 A. + 2tmA! + A" - 2* (imX -f X') + AS9 2) r cos r^ + XS9 2i , r sin r<£ = |-0 3t Wl , 

- m 2 X + 2(mX / -{- X" + 2* ((wA + A 7 ) + AS0 4| r sin r</> + XS9 6i r cos r^> = 0. 

Is a solution we now take 



V . (20) 



In these summations all the O's in the coefficients of p and or are to be included 
except O 1>0 and G 50 . A and X are constants, and the other coefficients functions 

of (/>. 

Now substitute these expansions for A and X in (20), and equate to zero the terms 
involving no except O 1>0 and 5u . We then have 



Whence 



-m 3 A -2«mX + e u A = : k®s,m> 
— m 2 X + 2/cr??^A u 4-O 5iU X = 0. 

Ao - |0 3 , m (O 5}0 -m 2 )™ {(G li0 -m 2 ) (e 5}U ^m 2 )-~4/cW} ? 
X = -~Kime 3m + {(Q h0 -m 2 ) (e 5i0 -m 2 )-4/cW}. 

Next, taking the coefficient of lr , we have the equations 

--m 3 A 1>r 4- 2imA / lfr + A // J)r --2«: (imX lir + X / 1)r ) + O 1)0 A l50 + A u cos r<f> — 
: ' -m 2 X l5r + 2imX / 1)r + X // 1)r + 2/c(^A 1?r + A / 1)r ) + e 5)0 X lir == 



>. (21) 



*' The use of m here to represent an integer is to be carefully distinguished from its previous use to 
represent mass. 
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We shall form the solution by taking only e tW> in the term cos r<p. Changing the 
sign of r will then give the other part. Assuming that A hr and X hr vary as e cr<j> , we 
have 

A ljr (-m 2 -2mr-r 2 -!-0 ])O )-2/ci(m + r)X ltr = — |A 0? 

A ltr {m + r) 2/a + X 1?r (-m 2 -2mr~r 2 + e M ) = 0. 
From these 



^ , . . . yZ/ZiJ 



and 



A 1>r = -|A {9 M - {m + rf} -S-[{e li0 - {m + r) a } {9 5 , - (m + r) 3 }-4«*(m + r) 8 ], 
X lir = + |A . 2,«(m+r)-h[{9 1 , -(m+r) 2 } {0 5 , o - (w + r) 2 }-4,c 2 (m+r) 2 ]. 



On determining the corresponding values for the term e - ""* and combining the two. 
we have 



A, r = -f A e"'* {9 M - (m + r) 2 } 4- [{e li0 - (m + r) 2 } {9 5>0 - (w + r) 2 } -4* 2 (m + r) 8 ] ) 



-iA e-^{e 3 ,o-(m-r) 2 }-=-[{e i , -(m-r) 2 }{e 5 , -(m-r) 2 }-4 K 2 fm-r) 2 ], 

r" (23) 
X,, r = A e"\i{m + r) -^[{O h0 - {m + r) 2 } {O bt0 - {m+rf } -i K 2 {m + rf] 

+ A u e- lH, Ki (m-r) + [{6 li0 - (m-r) 2 } {9 Bi0 - (m-r) 2 } -4k 2 (m-r) 3 ]. 

Expression (23) shows that A 1>r and X hr are factored by A , which is a multiple of 

9 3 , m . Now the terms in the expansions of A and X that we are seeking are A,,A, f 

and X ltr e lf r- Since both of these involve the product 63, ^i,,, it is clear that they 

may be neglected in comparison with the values of A and X . 

1 
We have further to determine the parts of A and X arising from a term — 9 6 , 1B e""* 

in the right-hand member of the second equations (4). These can be written down 
from the results already given, and are 

X = ^9 6 , m (9 5 , -m 2 )^{(9 1 , -7^ 2 )(9 5 , -m 2 )-4/cW}, 
A = «me 6 , m ^{(9 1 ,o-m 8 ) (9 6 , -w 2 )-4 K 2 m 2 }. 
Hence to the degree of accuracy we are using, we may summarise the results as : . 



2 [9 3jM! (9 5 , -m 2 ) cos m^ + 2rf»e 6>m cos ni^-f [(9 li0 -m 2 ) (9 5 , -m 2 )~4 K 2 m 2 ], 



m 



or 



= 2 [2 K m9 3 m sin m^ + 9 6 , m (9 5 , -m 2 ) sin m<p]^[{e h0 -m 2 ) (9 6 , -m 2 )-4 K W]. 



• (24) 



VI 



Except when the denominators are small, it is seen that, owing to the very small 
factors 9 3>m and 9 6 , m , the values of P and a- derived from the above equations are very 
small. 

YOL. COXXII.— A. 3 
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§ 4. Discussion of the Solutions of the 'Equations for the Case of Equal Particles. 

(i.) The complementary function. 

Equations (17) and (18) which determine the value of the exponent c, may be 
re-written here, 



. . (25) 



(26) 



-^©2 f 2» + ^Q 4f 2»— i" (ch,o-n 2 ) 5>2 J cos 2t 

2c(a lf0 e 5i0 -w 4 ) - {^(eo 5 o-^)Oi,2^^ 2 e 2}2ft 

+ ^ 2 4) 27l - |n (a lf - n 2 ) 5 , 2 J sin 2t J 

In these a lj0 is determined by the relation 

(a h0 -n 2 ) (e 5f o-V) - 4/cV. 

It is noticeable that the coefficient of sin 2r in (26) is n times that of cos 2r in (25). 
Owing to the smallness of the quantities (excepting 1O and 5fO ), it is clear that 
the coefficients of cos 2r and sin 2r are both very small quantities. Now real values 
of c are only given by real values of r, and conversely. Hence in order that (25) 
may give real values of r it is necessary that the expression 

{Oi,o-n 2 )(e^ ~-n 2 )-AK 2 n 2 (27) 

should be less than, or at most equal to, the coefficient of cos 2 T . That is, the real 
values of c will be in the vicinity of these values of K that make (27) vanish. The 
actual limits of the zone in which real values of c are found will be given by 

(0i,o™^ 2 )(O5,o"~^ 2 ) = 4KV±{i(0 5iO -n a )0 lt3 „ + icn0 3 ^ . (28) 

There are four groups of signs possible in this expression, and there will result four 
values of k. The outermost and innermost of these will define the zone in which 
some real value of c appears, and this zone will be the zone of instability. Owing, 
however, to the extreme smallness of the coefficient of sin 2r in (26), it is clear that c 
will be extremely small, in general; that is, the modulus of instability will be small 
and departure from the zone will be slow. In one case, however, c may be quite 
large. That is, when the coefficient of c, a 1?o 5)O -~n 4 , is exceedingly small. 

Each of the quantities is a function of ~~ f or of k. Further O h0 and 5>o involve 

a 

both the mass of the particles and the number of them. Both of these are entirely 
unknown. All that can be said is that Maxwell's criterion,* that is, 

v <~z, 

p 

* Tisserand, c Mec. Celeste/ vol. ii., p. 184. 
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where v is the ratio of the mass of a particle to the mass of Saturn and p is the 
number of particles in the ring, must be fulfilled. 

v appears in the expressions for 6 1}0 and 6 M in the form vL s . It has been 
mentioned that L, < 0'0l94n 3 for all values of s. Hence vL s < 0*0194yn 3 . For the 
present we shall regard vL s as a variable parameter and discuss the solutions relative 
to this parameter. 

In order to locate the zone of instability, we equate expression (27) to zero. 
Writing it in full, but omitting the term involving */,• which will be exceedingly 
small and will hardly affect the result, we find 

{{Z-vL s )K 2 + n 2 } {2vL/ + n 2 } -4k 2 u 2 = (29) 

This equation, regarded as involving an unknown quantity 71 2 Jk, is precisely the 
equation used by Maxwell to determine the condition of stability of the ring of 
particles when unperturbed by any satellite. The condition of the reality of u 2 /k 2 
leads to the upper limit for v just quoted. In our problem we may take the unknown 
quantity as K 2 fn 2 , and then assuming a value for i/L s , solve the equation. The values 
of k (for differing values of n) will give the position of the zones of instability of a 
ring of particles of mass and number assumed. Or, conversely, taking a position of 
instability, as shown by telescopic observations of the ring, we may determine the 
corresponding value of yL s , which establishes the order of value of the mass and 
number of particles at that point. 

I have found that the latter process leads to no satisfactory result, and hence I do 
not record the work. 

It is interesting to examine the meaning of the condition previously referred to, 
that the maximum instability is found when (^ 1G 9 50 — n 4 ) is approximately zero. On 
referring again to equation (28), it is clear that the broadest zone of instability 
will be found, owing to the extreme smallness of the last member, when 
(0i,o~~ ^ 2 ) (05,o~^ 2 ) —4tic 2 n 2 changes most slowly with k. This will occur when the 
equation (29) has equal roots. Equal roots appear when, by the variation of the 
parameter vL s9 k/u passes from real to imaginary values, or when # 

„L, = 0'039. 

This is the upper limit of the criterion previously quoted from Maxwell, and 
would imply that all the particles were of such mass and number as to be on the 
border-line of instability. 

When vL s has this value, we find that 

2vL s (3-vL s )K*-rf = 0; 
or 

e 5 , ei,o~^ - o. 

* Tisserand, loc, ciL, p. 183. 

s 2 
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Whence, by (ll), 

This right-hand member is of the same order of value as the factor of sin 2r in (26). 
In this case, then, c may assume a high value. But it is noticeable that only at the 
limit of Maxwell's relation is great instability to be found. 

When expression (27) has a value far from zero, either by virtue of the value of 
k/u or the value of vL s!) it is clear from (25) that 2t is imaginary and hence c is 
imaginary, the solution being stable. 

It might be inferred from this that if values of vL s were chosen such that 
Maxwell's relation were not fulfilled the effect of the satellite would be to stabilise 
what would otherwise be an unstable system. As pointed out already, however, the 
original equations and their solutions, as given here, simply give the motion of the 
particles in the vicinity of certain circles. In some cases the motion may be such 
that the particles depart rapidly from this zero circle ; this we have termed 
instability. In other cases the solutions may indicate that the particles never move 
far from the zero circle ; and this type of motion we have termed stable. But it is 
clear that if a small arbitrary displacement were given to each of the particles in the 
latter case, nothing in this paper precludes the possibility of their departure finally 
from the zero circle. That is, they may be again unstable. What we have found 
here is a series of orbits for the particles when subject to the attractions of Saturn, a 
satellite, and one another. Those in which the particles have large inequalities result 
in collisions with the neighbouring rings of particles and hence a complete departure 
from their former positions. Those which have no large inequalities and hence avoid 
collisions with neighbouring rings of particles may yet prove unstable when an 
arbitrary disturbance is further imposed upon them. 

(ii) The particular integral. 

In the expression (24) there appears a denominator of the form 

(e li0 -m 2 )(e 5i0 -m 2 )-^v (30) 

Here m takes all positive integral values including zero. When the conditions are 
such, therefore, that expression (30) is approximately zero, the term in the particular 
integral will become very great and departure from the orbit will result. This 
expression is the same as (29), which, it has been pointed out, gives the positions of 
the unstable solutions of the complementary function. It may therefore be said that 
all the unstable positions are in the vicinity of the zero values of (30), and the * 
following remarks apply equally to both parts of the solutions. 

Referring to the form (29) it is seen that there are two variables, k/u and vL 8 . For 
a given value of vL s there are in general two values of ic/n, and for a given value of k/u 
there are two values of vL s . In the figure (p. 125), the relation between K /n and vL 8 
is shown graphically, only those values of vL s which satisfy Maxwell's criterion being 
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chosen. It will be seen that k/u increases slowly from unity as vL s increases from 
zero, until vL s reaches the value 0*039. At this point the curve turns back and rises 
rapidly to an asymptote at vL s = 0. 

In the case when expression (80) is exactly zero, it is seen from elementary principles 
that the independent variable <f> would appear explicitly. With passage of time, 
therefore, p and a- would increase linearly in magnitude and there would be complete 
departure of the particles from the vicinity of r = a. 

§ 5. Case where the Particles forming the Ring are of Unequal Masses. 



The previous equations (2) were reduced to the form (3) on the supposition that all 
the masses m k were of the same value m, = vM.. We now proceed to the modifica- 
tions introduced when these masses are all distinct in value. 

Equations (2) with the same reductions as before, but maintaining the separate 
values Wx, become : 



p A — Zk<T k 



/I Ait] 



"\ 



/c YV \ -j- ( _A + . . . + hi cos id> + ...)- cos <p \ — /c 2 E ; 

[da 



+ 



d 2 

Sic 2 + A' V y— (%b + . . . + bi COS id> + . . . ) — ac 2 F ; 

da J 



-_I 



Pa + 2/c Gr^x/V 



fi 



-f [ A /4 V ~t~ (&i sin <p + ... + ibi sin i<p •+• . . . ) — A /4 V sin ^-/c 2 HJ o- A 



H 



<r" x + 2k p\ = AY 2/ * (& L sin c/> + . . . + ib { sin ^ + ...) — AY 4/s sin -f j^E 7 , 



f 



/ 2 /4/ s 



i/WM ... -m 



• 



(2 



a 



sin i^ + . . . j + vkk*!* ( . . . + i&i sin i<p...) 



/c r x 



<Tx+2/C*G /Xj x/J /x 



/x 



+ [ — AY 2/s ( . . . + i% cos ^ + . . . ) + AY 2 cos - /c 2 IT A ] o-x 



+ 2/C J u.XCT.u» 



/x, A^ /x< 



/x 



> 



In these equations 



Ex = 2 ^ 



n 



M [4 sin (m — A) tt/ 



r^ 



+ cos (/x— A) 2irjn\> 



■ (31) 



R = Z m * 



1 



3 



/* 



M [8 sin 3 (/x— X) 7r/r^ sin (/x— A) 7r/?i. 



G 



/.<,. x 



3 

8 



m /x ( cos Q*— \) 2tt/^ ___________ 

M [8 sin 3 (}x — \) 7r/n sin (/x~ A) irjn 



+ 2 cos (a—A) 2irjn 
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H, 



fX 



f COs( M -~\)7r/n . / x 9 , 

,, i" — r-p 7 r 7 - + Sill (/x — A) Ztt/71 



/_, — — — 



J 



ix, A 



m /x [ cos (,u—\) 7v jii 



M 18 sin^/x-A)^ 



+ sin (/ul — \) 2ir/n 



r 



Ef 



F, 



y m^ ( cos (m-X) ttM . , x 9 / 1 

Zj -^ S- < ~-~^—~~~—~~— — sin ( /x — A ) Air hi r •> 

t, M (4sm a (/x-A) 7r/n A ! ' 



m^ [ 3 cos (/x^X) 7r/W 
"M"L8sin 2 (^-x)7r/^[' 



jX 



G f 



m 



/.i, A 



^ 



M 



— ~^4t \nn + 2 sin (/x-x) 2rr/n\ , 

8 sm 2 (juL—\)>7r/n J 



tt/ - t m ^ / cos (^~ A ) 2?r M A sin 2 (/x-~X) 2tt/u / x 9 / 

XT A ~~ -^.~i\ar 1 -^-T^'TY~^^~'~^ m Y^^^ "~~ """ """** =^r"nTT™™~^^ ~"~ COS ^u — A / ZiTT/ tl 






' M .18 sin 3 (yu — X) ?r/n sin 5 (ju — X) wfn 



3 f 



ix, A 



■ m ^ [ cos (/u— -\) 2?r/^ __ -^ sin 2 (/x—X) 2ir/n 
M 1.8 sin 3 (/x — X) w/r& sin 5 (/x~X) ?r/^ 



cos (jul — X) 2irjn\. 



9 * 



(32) 



These equations may be written, with a slightly different meaning of the quantities 
from those formerly obtaining, in the form, 



p\ — 2kot\ + p A 2 {0 1>r cos r<j>} + /c 2 2G^ xp/x + o-a2 {0 2 , r sin r^} 



■\ 



/x 



7' 



+ /c 2 2 J^ k a-fx = 20 3 , r cos r 0. 



^ 



o- /; A ■+• 2/c// A + p A 2 { 4j v sin r^ } + r2G^ t A p /X + o^S { 5f r cos r<^> } 



y. (33) 



/x 



•i- /c 2 S J / /x> A o- M = 20 6ir sin r^. 



IX 



r 



y 



We shall determine the particular integral arising from one term of the right-hand 
member of the first equation, writing it typically ^0 3> m e im4> . 
Assume that 

Pa — -A-a^ 5 



and 

for all values of X. 

Equations (33) then become 



<r x = X x e tm</> 



■ m 3 A A + 2miA\ + A 



// 



- A -r^«/&xi. A -r^i. A™2/c(/wX r }-X / A ) + A A 2 {0 ls7 , cos r^ 



+ /c 2 A^Cx^ A + X A 2 \02, r sm T<j>\ + /c^S-X.^ J ^ A = 2*0,3 m , 



■m 2 X A -f 2imX\ + X /; A + 2/c (mA A + A 7 A ) + A A 2 {0 4) r sin r^} 

-f-^SA^G^x-j-XxS {0 5 , r cos r^} 'f/SX^J^x = 0. 



. (34) 



There are n pairs of equations in this form corresponding to the n values of X. 
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As a solution we now take 

A • .XV.Q ~t" 4~-l<£*uJ-±~tf g\Is y $ t • » • j 

with the same restrictions as before. 

Substitute in (34) and equate to zero the terms involving no 9 but 1>o x and 5tO A . 
Then, for all values of A from 1 to n, 



lib _OLa 



2 K tmX$ +e x h0 A£ +« 2 2;AjG^x + ic a SXf J„ iX 



la ^ 

2 ^3, m) 



^< 



ja 



m 2 xj +2«wAi +85,0 X£ +/c 2 s;At; gv,x+k 2 sXs- J'^x = 0. 



> . . (35) 



[X 



fX 



These 2n equations can be solved by the usual processes to give the values of the 
constants Aq and XJ. It is not necessary for us to work out the results in detail, it 
is sufficient to note that the determinant of the left-hand members will appear as the 
denominator in each case. The determinant is the following : 



2/am ' , 

2fif 

K x: * 1, 2 •> 
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• (36) 



J 



This determinant corresponds to the denominators in expressions (24). When it 
vanishes or becomes small, it is clear, as before, that the terms of the solution tend to 
become large, and instability follows. 

In estimating the values of F, G, H and J, which appear in the above determinant, 
it is to be noted that mjM. is exceedingly small for all values of //. But the quantities 
in which it appears may be large by virtue of the small denominators which are 
involved. In the expression for F x , the term 3/sin (/* — A) irfn may be neglected in 



comparison with the first term for large values of n. Also, X 

1 



, m 



a 



will 



m 



M 8 sin 8 (n — A) 7r/n 

lie between zero and ~ 2 — -7-—-^ — ■ T -— T - since all the signs are positive, if m is the 

M 8 sitr* (m-A) irjn to l 

greatest value of m lx appearing in the ring. Hence F A lies between zero and 
0*0096 w 3 m/M in value. 

In the same way the value of G^ )A will arise almost wholly from the first term. 
The largest value it may have will be mri A f '8tt 3 M or 0'004n ;J, m/M . E A , H X5 J^a, E^, F\ and 
G f ^ x are seen to be one order lower in the reciprocal of sin (pc— A) tt/u and therefore 
may be neglected. H ; A has the limit — 0*01927?/ 3 m/M, and J ; /X( A the limit — 0'008n 3 m/M. 
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We shall assume that the number of particles in any ring is large. It is probable 
that they vary in magnitude from the infinitesimally small up to the limit given by 
Maxwell. Hence the values of the expressions F A , G^ A , H\ and J\ will vary over 
a range of values, between the given limits, as X takes its successive values. 

Reverting to determinant (36), we see that it may now be written 
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(37) 



7 



For all conditions satisfying Maxwell's criterion, the quantities G x , /X3 J ; x,^ will be 
small. So that, provided k 2 is not too great, the value of the determinant (37) will 
be small for those values of K that satisfy the relation 
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(38) 



This relation is satisfied by those values of k which satisfy the equation 



(-m a +e\ ){-m*+e\ )-4 K 



9 9 

m 



U, . . . , . , [Otj J 



where X takes all its integral values in turn. Further it is easily shown that, on any 
distribution with n large, -F A = ^H ; A . Hence we fall back upon the same type 
of equation as we had in the case of equal particles (equation (29)) where we replace 
vh s by F A . 

Instead of treating the equation (39) separately for the various integral values 
of X, sinre n is large, we may imagine a single equation with the assumption that 
F A is an arbitrary variable parameter. The determinant (37) will then be small, and 
instability result for all the values of K given by (39), for all values of the parameter 
F A that exist. With a wide range of values of F A corresponding to a wide range in 
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the magnitudes of the masses of the particles, we may expect to find a broad region 
of instability. 

It can readily be shown that the condition (39) would also be produced if the 
general case of unequal particles were solved for the complementary function in the 
same way as has been done for the case of equal particles, which produced (29). 
This work is not reproduced owing to the length and complexity of the expressions, 
and also because the results are wholly contained in the condition (39) produced from 
the particular integral. 



6. Application of the Results to the Saturnian System. 



Equation (39) written out in full is 



{/c 2 (3-F A ) + m 2 } {2F A K 2 + m 2 }-4/c 2 m 2 = 0. 



(40) 



In this equation m is any integer and F A may vary between zero and 0'0096n 3 m/M. 
As the distinctions indicated by the suffix A are 
now of no importance, it may be dropped. The 
solutions of (40) will give approximately the 
positions where divisions in the Ring of Saturn 
may be expected. 

For any given value of F, there are four 
values of tc/m, two pairs equal with opposite signs. 
For any given value of tc/m there are two values 
of F ; one, however, being greater than the 
Maxwell limit, is excluded. The limiting value 
of F for real values of /c/m is 0'039. This is, 
of course, the same result as that found by 
Maxwell. 

The relation between /c/m and F is shown in 
the figure, and the table shows actual numerical 
values. 

We may readily assume that in the existing 
Rings of Saturn there are particles of all masses 
from the infinitesimal to Maxwell's upper limit. 
These will give rise to varying values of F, 
depending upon the masses of the particles 
adjacent to the particle under consideration. The 
maximum value of F is itself small compared with unity ; we shall then arrive afc a 
limit of k by taking F = in equation (40). We find thus that the boundary of a 
division should occur at k/tu = 1, for each integral value of m. 
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= 1 
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a\a! = 0, 
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= 2 


I, a)/&> = 2 , 


a/a 7 = 0*62996, 






m 


= 3 


j 60/CO = 2" , 


a/a 7 = 076289, 






m ■ 


= 4 


j *»/«' = | , 


a/a 7 = 0*82524. 





[It should be remarked that a positive value of /c gives positions without the 
satellite orbit, and a negative value of k gives positions within. As k appears in (40) 
in the form of a square, both positive and negative forms result. We should 
therefore have the same phenomena in a ring of particles beyond the satellite orbit 
as we find within]. 

The result a = implies a division of the ring at the origin. This would fall 
within the planet itself. But if the zone consequent upon the variation of F is 
extensive, it may extend beyond the surface of the planet and show a clearance 
of particles there. 

For m = 2, a\a r = 0*62996. In the case of satellite Mimas this should indicate 
the commencement of a division in the ring at distance 16'9 ;/ . Cassini's Division 
begins at 16 '87" and ends at 17'64 // . This agreement is very remarkable. 

Reference to the figure shows that in the vicinity of ic/m = 1, k increases with F. 
But as k increases so does a/a'. Hence the instability caused by the larger values 
of F should be in positions corresponding to larger values of ^ that is, to larger values 
of a/a\ In other words, the division should extend outwards. This agrees with the 
observational data just quoted. We may then attribute the production of Cassini's 
Division to Mimas. 

For n = 3, a/a f = 0*76289. 

For satellite Mimas, this should cause a division at distance 20'46 // . This is just 
beyond the outer edge of Ring A, which terminates at 20*0 1 ;/ . 

Considering next the satellite Enceladus, we should find a division at the origin 
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for n = 1, and at distance 21*69" for n = 2. The last is again just beyond the limits 
of Ring A. 

The remaining satellites all produce instability at the origin, but the other points 
at which this occurs are outside the existing ring. 

We may use the observations of the dimensions of Cassini's Division to determine 
the maximum value of F appearing. As we have already found, the inner edge 
corresponds closely to F = 0. The radius of the outer edge is a = 17'64". Hence 
for satellite Mimas ajaf = 0*65781, giving KJm = 1*0720. 

If now equation (40) be solved for F, taking this value for it/m, the result is 
F = 0*0173. Hence we may conclude that F ranges from zero to 0*0173. 

Using this value of F, we now proceed to the study of the roots of equation (40). 

Solving, we find 

K /m = ± 1'0720 and K /m = ± 2*8917. 

Take m = 1. Then 

a/a' = 0*1712 and a\a! = 07535. 

We may expect to find a clearance of particles from a\a ! = to 0*1712; and from 
aja f = 0*7535 to unity. 

The first gives the extent of the clearance near the origin. 
For the various satellites its dimensions are : 

Mimas a = 4*59", 

Enceladus a = 6 '16", 

Tethys ......... a = 7*30", 

Dione a = 9*34", 

Rhea a = 13*07", 

Titan . . a = 29*94". 

The radius a = 9*34" indicates approximately the inner radius of the Cr6pe Ring, 
while a = 13*07" indicates more closely the inner radius of Ring B. 

Applying the second ratio, aja' = 0*7535, to Mimas, we find radius a = 20*2" 
There should be a clearance of particles from 20*2" up to the satellite itself. This 
indicates with considerable precision the termination of the whole ring, which has a 
radius 20*01". 

These results are subject to modification owing to the effect of the oblateness of 
the planet Saturn and the influence of one ring upon another. But the agreement 
of theory and observation in this first approximation is sufficiently remarkable. 

The interpretation of the effect of Dione and Rhea on the inner parts of the ring 
is not clear. From the theory one would expect that any satellite could affect a 
clearance of particles from the origin up to a radius given by a/a f = 0*1712. In that 
case Titan, the largest of the satellites, should dissipate the whole of the existing 
rings, for this ratio carries us far beyond the outer radius. 

T 2 
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There are therefore two facts to explain. First, the existence of the Cr6pe Ring 
within the dissipative area of Rhea, and second, the existence of the bright rings 
within the dissipative area of Titan. In connection with the first, Lowell has noted 
a definite black band within Ring B, so that there is a clearance of particles between 
the Crfrpe Ring and the bright rings. It would appear as though the dissipative 
power of the satellites was only effective near the outer boundary of the unstable 
area about the origin. To discuss this, let as examine the analytical results. 

It has already been pointed out how very small the exponent c is, as given by (18), 
indicating a very slow rate of dispersion. Consider, instead, the numerators of the 
expressions (24), the vanishing of the denominators of which causes the instability. 
The numerators are small because of the quantities 9 3>m and 6?m . In the case under 
discussion, m = 1. From (5) 

Ida 
Using the well-known expression for 6 l5 # we find 

e 3 , 1 = /A' 1 Mt* 2 +!f-t« 4 +...} 
= ,/{»7«-i}- 8 {K+M* 6 +-.-} 

9 6il = A 2 {a{a + ia s + T %^a\..)-a 2 } 

= /{ w y w -i}- 2 {f« 4 +T¥%« 6 .-.}. 

For small values of a, w'/w is small, and the value of (&//«— 1)~ 2 will be greater than, 
but not far from, unity. Hence the values of 9 3>1 and 6>1 depend approximately 
upon the fourth power of a or a/a'. It is clear then that the numerators in (24) will 
be vanishinffly small except for the larger values of a/a'. 

The' physical meaning is that, while instability will always take place when the 
denominators vanish, the rate of dissipation will be small except for the largest values 
of a which are permissible. There will also be a uniform grading in the rate of 
dissipation as a increases. 

Applying this result to the case of Saturn's satellites, we may expect to find 
actually a clearance only near the outer limits of the areas under consideration. The 
areas of clearance of the first three satellites fall within the body of the planet. 
Dione causes the clearance between the surface of the planet at 8'65 /; and 9*34", 
which is approximately the commencement of the Cr£pe Ring. The limit of the area 
of clearance of Rhea is 13 '07", and only near that boundary is the action effective, 
the Crepe Ring being undispersed in the weaker part of the field. The bright rings 
are clearly in the weak part of Titan's field of clearance, and so continue to exist. 
It is obvious, however, that with passage of time the Crepe Ring will be dispersed 
by Rhea and the whole by Titan. 

* Tissehand, vol. i., p. 272. 



SATELLITES UPON THE FORM OF SATURN'S RING. 129 

We have found the maximum value of F appearing as 0*0173. It was previously 
shown that the limiting value of F was 0'0096^ 3 m/M. Hence 

0'0096^ 3 m/M = 0*0173, 
which gives 

- /M 0-0173 



0'0096^ 3 
= 1'8/^ 3 . 



That is, the size of the largest particles is just below that given by Maxwell's 
criterion. 

§ 7. Summary and Conclusion. 

( 1 ) Assuming that a planet is surrounded by concentric rings of particles performing 
approximately circular orbits when unperturbed, and that the influence of one ring 
upon another may be neglected to a first approximation, the effect upon these rings 
of a satellite performing also an unperturbed circular, orbit is discussed. 

If the particles in the rings are all equal, it is shown that we should expect, 
in certain places, large perturbations to take place, such that the particles in a 
particular ring would leave that ring and mingle with those of other rings, and so 
leave a "division." 

(2) As there is no reason to believe that the particles in any ring are all equal, the 
analysis is extended to cover the case of unequal particles. 

We assume that in any ring the number of particles is large, and that therefore we 
shall probably have one specimen at least of all particles from the smallest to 
the largest. 

It is then shown that the divisions would become more extended, and therefore 
more readily visible. 

(3) On the supposition that some of the particles at any rate are indefinitely small, 
we obtain Cassini's Division at once. On making use of the dimensions of this 
division to estimate the greatest magnitude of the particles in any ring, we find the 
following results : — 

.Satellite Mimas should produce a clearance of particles from radius 20*2 // up to 
itself. The ring should therefore terminate at 20'2 // . Observation shows that it 
terminates at 20'01 // . 

Satellite Mimas should produce a division from radius 16'9 /; to 17'64 /; . (This last 
measurement was used as a datum for estimating the magnitude of the greatest 
particles.) Observation gives the limits of Cassini's Division as 16'87" and 17*64 // . 

Satellite Dione should produce a clearance of particles from the region of the 
surface of the planet up to radius 9 '34". The Cr6pe King is observed to begin with 
a diffused edge at 10* 8 3". 
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Satellite Rhea should also produce a clearance of particles up to radius 13*07 . 
The inner edge of Ring B is observed to commence at 13 '2 1". 

The existence of the Crepe Ring in a dissipative area is also discussed. 

(4) By the inclusion of the effect of the oblateness of Saturn and the influence of 
one ring of particles upon another these results might be still further improved. 

The theory presented therefore gives a closely quantitative account of the salient 
features of Saturn s Ring. The numerous smaller divisions observed by Lowell and 
others are not accounted for ; but, for the reasons given in § 4, their existence is not 
excluded. 

(5) The dimensions of Cassini's Division show that particles of all sizes up to a 
limit just short of that imposed by Maxwell for stability exist in the rings. 



Appendix on the Data of the Problem. 

I. Dimensions and divisions of the ring in seconds of arc at mean distance^ 

Distance from centre of planet to — 
Inner edge of Cr^pe Ring . 
Inner edge of Ring B . 

Bl . 



Divisions of Ring B -i 



B2 
B3 
B4 
B5 
B6 



Outer edge of Ring B . 

Inner edge of Ring A . 

Division in Ring A . . 

Outer edge of Ring A . 

II. Equatorial diameter of Saturn . 



III. Elements of satellites ;• 



Mimas . 
Enceladus 
Tethys 
Dione . 
Rhea . 

X 1 u<3jll < . 



Mean distance. 

26 '82" 
34-43" 
42-66" 
54-59" 
76-38" 
174-8" 



10-83" 
13-00" 
13'39" 
14-04" 
1474" 
15'32" 
15'69" 
15-95" 
16-87" 
17-64" 
19-00" 
20-01" 



// 



17-30 



Mass 
as fraction of Saturn. 

7. 10~ 8 

25. 10" 8 

11 . io- 7 

18 - 7. 10- 7 

4. IO" 6 

2-1 . 10-" 



* Lowell, 'Observatory Bulletin,' No. 68, and "Lecture" on April 26, 1916, in 'Journal of Eoyal 
Astron. Soc. of Canada.' 



